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Light element synthesis in supernovae through neutrino-nucleus interactions, i.e., the i^-process, 
is affected by neutrino oscillations in the supernova environment. There is a resonance of 13-mixing 
in the 0/C layer, which increases the rates of charged-current z^-process reactions in the outer He- 
rich layer. The yields of ^Li and ^^B increase by about a factor of 1.9 and 1.3, respectively, for a 
normal mass hierarchy and an adiabatic 13-mixing resonance, compared to those without neutrino 
oscillations. In the case of an inverted mass hierarchy and a non-adiabatic 13-mixing resonance, 
the increase in the ^Li and ^^B yields is much smaller. Observations of the '^Li/'^^B ratio in stars 
showing signs of supernova enrichment could thus provide a unique test of neutrino oscillations and 
constrain their parameters and the mass hierarchy. 

PACS numbers; 26.30.-|-k, 14.60.Pq, 25.30.Pt, 97.60.Bw 



A tremendous number of neutrinos are released from 
a core-collapse supernova (SN). These neutrinos inter- 
act with nuclei in the surrounding stellar envelope and 
thereby affect the synthesis of new elements. This so- 
called I'-process may be a major contributor to the pro- 
duction of several hght isotopes, such as ''Li, ^^B, -^^F, 
as well as a few heavy isotopes, such as ^^^La and ^*°Ta 
[1 H H H H EQ' However, the yields of these iso- 
topes may depend on the effects of neutrino oscillations, 
which was not taken into consideration in the above cited 
studies. 

Recent neutrino experiments on atmospheric 01, solar 
, and reactor neutrinos 0, significantly constrain 
most of the neutrino oscillation parameters. However, 
only an upper limit on ^13 is obtained |lCll | and the mass 
hierarchy remains unknown. Theoretical studies of neu- 
trino oscillations in SNe have been used to suggest po- 
tential constraints on ^13 and the mass hierarchy based 
on observed SN neutrino spectra. These studies indicate 
that the neutrino spectra from SNe strongly depend on 
6*13 and the assumed mass hierarchy [13. fisj . When the 
resonance of the 13-mixing is adiabatic, substantial con- 
version i>c ^ v^,T occurs in the 0/C layer for a normal 
mass hierarchy and conversion Vc v^^t occurs for an 
inverted hierarchy. These direct methods are of course 
limited by the fact that nearby core-collapse SNe occur 
rarely in the small detection volume given by current de- 
tector sizes and methods. 

Here we suggest an alternative method to study the ef- 
fects of neutrino oscillations, by considering light element 
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synthesis in SNe. Neutrino energy spectra change as they 
are transported through SN ejecta [IJ. This change will 
affect the production of light elements via the i/-process. 
The thermal neutrinos emitted from a cooling protoneu- 
tron star have a well-known, but not yet rigorously estab- 
lished, energy hierarchy; {e^^) < (spj < (ei/^,,,p^,^) (e.g., 
[r3|'). Neutrino oscillations could thus increase the aver- 
age energies of Vc and D^, and consequently the rates of 
charged-current z/-process reactions could be much larger 
than expected from models without oscillations. There- 
fore, the yields of the light elements may increase signif- 
icantly. 

We investigate nucleosynthesis of light elements ^Li 
and "'^^B through the ^/-process in SNe taking neutrino 
oscillations into account. Since the other j/-process el- 
ements are mainly produced in the 0-rich layers 0, , 
they are not expected to be affected by neutrino oscilla- 
tions. The ^Li and ^^B yields in SNe can thus be used 
as probes of neutrino oscillations. We show that the ^Li 
yield increases significantly through neutrino oscillations. 
The dependence of the ^Li/^^B ratio on the mixing pa- 
rameter, 013, provides an observable signature that could 
be used to constrain its absolute value and the neutrino 
mass hierarchy. 

Neutrino luminosities are assumed to decrease expo- 
nentially in time, with a decay time scale of '--^ 3 s 
0, 0, 0, S IE ■ The total neutrino energy is assumed 
to be fixed at 3 x 10^^ erg. The neutrino energy spec- 
tra at the neutrino sphere are approximated with Fermi- 
Dirac (FD) distributions with zero chemical potential. 
The neutrino temperatures of i^c, v^, and {v^i^t and 
are set to be 3.2, 5.0, and 6.0 MeV as adopted in j^. 
These energy spectra change during the subsequent pas- 
sage through the outer stellar layers by neutrino oscilla- 
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tions. 

In order to evaluate the reaction rates of the ^^-process, 
we need the cross sections as functions of neutrino energy 
because the spectra changed by neutrino oscillations no 
longer follow the FD shape. We assume that the cross 
sections of the charged-current reactions of "^He and ^^C, 
He(i/e,e-p)3He, '^}ie{De,e+n)^R, i2c(j/e, e-p)"C, 
C(i>e,e+n)"B, i2c(z.„e-7)i2N, and ^^C{P,,e+jy^B, 
are expressed as a power law a{e,y) ~ o-o(£i/^eth)", where 
eth is the threshold energy. Coefficients of the functions 
are determined such that the reaction rates deduced us- 
ing these cross sections (and assuming FD energy distri- 
butions) fit the rates tabulated in [I^. Details are pro- 
vided in 0. For the other i/-process reactions, we use 
the reaction rates with FD distribution of the neutrino 
spectra. 

Recent neutrino experiments |^ 0, 0, have de- 
termined most of the values of the mass squared dif- 
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and sin^ 26*12 = 0.816, 
< 1 X 10"^ The posi- 
tive value of Ato|i corresponds to "normal hierarchy" , 
i.e., TOi < TO2 < TO3 and the negative value corresponds 
to "inverted hierarchy", i.e., TO3 < TOi < TO2. We nu- 
merically solve the mixing probabilities of neutrinos for 
each neutrino energy by Runge-Kutta methods and using 
the exact solutions of the oscillations described in [l7l |. 
By convolving the mixing probabilities and the neutrino 
spectra at the neutrino sphere, we evaluate the neutrino 
energy spectra taking neutrino oscillations into account. 
We do not include CP phase S. Based on |0| we assume 
that the effect of CP violation will not be seen because 
the spectra of and Vr {i^t) emitted from the neu- 

trino sphere are the same. The change of the spectra due 
to oscillations is calculated using the density profile of a 
presupernova model. 

We use the same SN explosion model as in Q ■ The 
presupernova model is the 14E1 model constructed for 
SN 1987A in [l^. The SN explosion is calculated using 
piecewise parabolic method code |23, |^ ■ The explosion 
energy and the location of the mass cut are set to be 
1 X 10^^ erg and 1.61 Mq. The detailed nucleosynthesis 
in the SN is calculated using a nuclear reaction network 
including 291 species 0. 

Figure 1 shows the mass fraction distributions of ^Li 
and "'^'^B in the SN ejecta with neutrino oscillations of 
sin^ 2013 2 X 10"^ and for those without oscillations. 
In the case of a normal hierarchy, the mass fraction of 
^Be with the neutrino oscillations is much larger than 
that without oscillations in the He layer. There is a 13- 
mixing resonance for neutrinos in the 0/C layer and the 
resonance is adiabatic in this case. Thus, the energy spec- 
trum of i/c in the He/C layer becomes almost the same 
as that of Vfj^^r in the 0-rich layer. Beryllium 7 is pro- 
duced through ^He(j/, i^'n)^IIe(Q;, 7)^Be. Owing to the 
neutrino oscillations, the reaction rate of *IIe(i'e, e-p)'^IIe 
becomes larger than that of *He(i^, z/'n)'^IIe. The 



10 



J8 10-9 



10 



10 



O-rich 



0/C de/C He/N 



^1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

7ia-> /^^- JLi 
Be/ ^1?=:$^ 




r 


f 4 


: 








































FIG. 1: The mass fraction distributions of ^Li and its isobar 
^Be (upper panel), and ^^B and ^^C (lower panel) in the case 
of sin^ 2^13 = 2 X 10"^. Thick lines indicate the distributions 
of ^Li and ^^B. Thin lines indicate the distributions of ^Be 
and ^'^C. Solid lines and dashed lines correspond to a normal 
hierarchy and inverted hierarchy, respectively. Dotted lines 
correspond to the case without neutrino oscillations. The 
horizontal axis is the interior mass in units of the solar mass. 



mass fraction of ^Li including neutrino oscillations is 
also larger, but the increment is much smaller than 
that for '^Be. The main production process of ^Li is 
^116(1^, i^'p)"^H(q;, 7)^Li and the corresponding charged- 
current reaction is '*IIe(i/c, e~'"n)'^H. However, there are 
no resonances for antineutrinos. 

The effect of neutrino oscillations on the mass frac- 
tion distributions of ^^B and ^^C is similar to that 
for ^Li and ^Be. The mass fraction of ^^C with the 
neutrino oscillations is larger than that without oscil- 
lations in the He layer. During the i/-process, ^^C 
is produced through ^^G{iy, i^'n^^C and partly through 
i2C(^c,e-p)"C. The reaction rate of ^20(1/^, e-]5)"C 
with the oscillations becomes larger than that without 
oscillations by about one order of magnitude. The mass 
fraction of "'^^B with oscillations is only slightly larger 
than that without oscillations. The main production pro- 
cess of "B is ''He(z^,j/'p)^H(a,7)^Li(a,7)"B. The cor- 
responding charged-current reaction is *He(t'o, e^n)"^H. 
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About 12% - 16% of in the He layer is produced 
from through ^^C{i^,iy'py^B and i2c(p„, e+n)"B. 
The ^^B abundant region in the He layer is inside the ^Li 
abundant region because of the decrease in peak shock 
temperature as one moves outward in the star. The in- 
crease in the production through ^^C{Dc,e^n)^^B is 
not as large because of the absence of resonances for an- 
tineutrinos, as mentioned above. In the 0-rich layers, 
light element production is not influenced by the neutrino 
oscillations. The oscillation amplitude in these layers is 
too small because of high densities (e.g., 

In the case of an inverted hierarchy, mass fractions 
of ^Li and ^^B are larger than those for a normal 
hierarchy. The reaction rates of "'He(f'c5 e+n)'^H and 
^^C{Dc, e+n)iiB become larger owing to an adiabatic res- 
onance of ^fj.,T- However, the increment of the mass 
fractions of ^Li and ^^B is less pronounced than that of 
^Be and ^^C for a normal hierarchy. This is because 
the average energy of Vc is larger than Vc at the neu- 
trino sphere and, therefore the difference from the aver- 
age T- {Vfj^^r) energy is smaller. On the other hand, the 
mass fractions of ^Be and ^^C are slightly larger in the 
mass range Mj. > A.^Mq and slightly smaller inside the 
range of the He layer. There is no 13-mixing resonance 
for i/e, so that substantial conversion of ^ v^.t does 
not occur. At the same time, some ^Be and ^^C capture 
neutrons produced through *He(f'e, e^n)'^H. 

Figure 2 shows the ratios of the ^Li and ^^B yields 
with neutrino oscillations in comparison to those with- 
out oscillations, hereafter called yield ratios, as a func- 
tion of sm^29i'i. The yields of ^Li and ^^B without 
the oscillations are 2.36 x IQ-'^Mq and 6.26 x IQ''^ Mq. 
The yield ratio of ^Li is at most 1.88 in the case of 
sin^ 2013 > 2 X 10^"^ and normal hierarchy. This in- 
crease in the yield is due to the adiabatic 13-mixing res- 
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FIG. 2: The yield ratios of ^Li and ^^B with the relation of 
sin'^ 26113 • Thick solid line and thick dotted line are the yield 
ratio of ^hi in the cases of a normal hierarchy and inverted 
hierarchy, respectively. The thin solid line and thin dotted line 
are that of ^^B in the cases of normal and inverted hierarchy. 
The case of sin^ 26ri = is also calculated (see text). 



onance. In the case of 2 x 10~^ < sin^ 26113 < 2 x 10~^, 
the yield ratio of ''Li increases with sin^ 'iOi^ . In this 
013 range, the 13-mixing resonance changes from non- 
adiabatic to adiabatic with increasing in sin^ 20i3 . In 
the case of sin^ 26*13 < 2 x 10~^, corresponding to nona- 
diabatic resonance, the yield of ^Li is about 1.13. In the 
case of an inverted hierarchy, the dependence of the ^Li 
yield on sin^ 20i3 is similar to the normal hierarchy case. 
However, the increment of the yield ratio is much smaller. 
The smaller difference of the average energy between Vf. 
and i>p_^T reflects the smaller increase in the yield ratio. 

The dependence of the ^^B yield ratio on sin^ 20i3 is 
similar to that of ''Li. The ^^B yield ratio is about 1.27 
even in the case of adiabatic range of ^13 and normal hier- 
archy. This value is much smaller than that of ^Li. Neu- 
trino oscillations raise the rate of ^^G{vcTe~pY^C and 
the yield. However, the increased ^^C yield is still 
small for the total yield of "'^'^B. In the inverted hierarchy 
case, the maximum yield ratio of ^^B is slightly smaller 
than that in the normal hierarchy. In this case the contri- 
bution of ^^G{i'e,e^nY^B and '*He(i^o, e'''n)'^H increases. 
As shown in ^Li case, however, the increment is small 
due to a small difference of the average energy between 
Vc and v^^T- In the limit of sin^ 26*13 = 0, the ^Li and ^^B 
yields still slightly increase due to the residual mixing 
other than 13-mixing as shown in Fig. 2. 

We solve for neutrino energy spectra changed by neu- 
trino oscillations in the density profile of a presuper- 
nova star. We expect that the influence on the spectral 
changes due to neutrino oscillations caused by the pass- 
ing shock is small. When the shock front is in the 0-rich 
layers, the density behind the shock front is still so high 
that the shock wave does not affect the oscillations. Af- 
ter the shock front passes through the 0/C layer, the 
change of the density profile affects the mixing proba- 
bility of neutrino oscillations. However, most neutrinos 
have passed before the shock arrival at the 0/C layer. 
Details are discussed in [T^ . 

In our previous studies 0, 0; we constrained the 
neutrino temperature with Galactic chemical evolution 
(GCE) arguments (e.g.,j22|) for ^^B. However, the re- 
maining model uncertainties still render the effects on 
the observed ^Li and ^^B abundance trends from neu- 
trino oscillations somewhat ambiguous. Still, the possi- 
bility to obtain, or at least constrain, fundamental neu- 
trino properties from these observations encourage us to 
pursue these arguments further. We consider the depen- 
dence of the ^Li/i^B ratio on sin^ 20i3 taking account un- 
certainties of neutrino energy spectra. We consider two 
additional spectral parameter sets: {T^^,Tp^,Ti, ^,E^) 
= (3.2 MeV, 5 MeV, 6.6 MeV, 2.4 x lO^^ gj-g) ^nd (3.2 
MeV, 4.3 MeV, 5.2 MeV, 3.5 x 10^^ ^^.^y ^he obtained 
i^B yields for these two cases without neutrino oscilla- 
tions are 7.3 x 1O~^M0 and 3.3 x 10~'' Mq, correspond- 
ing to the maximum and minimum values satisfying the 
GCE models for ^^B 0. The corresponding ^Li yields 
are 2.9 x 10"'^ Mg and 1.3 x IQ-^M©. 

Figure 3 shows the number ratio of ^Li/^^B as a tunc- 
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FIG. 3: The number ratio of ^Li/^^B with the relation of 
sin^ 2^13 . The shaded ranges include the uncertainties of neu- 
trino energy spectra deduced from the calculations using three 
sets of neutrino temperatures and total neutrino energies (see 
text). 

tion of sin^ 2^13. The uncertainty due to neutrino spectra 
is included as shaded regions. The ^Li/^^B ratio in the 
case of adiabatic 13-mixing resonance and normal hier- 
archy is larger than that without neutrino oscillations, 
even with the spectral uncertainties included. Thus, the 
enhancement of observed ^Li/^^B ratio may constrain 
the lowest value of ^13 and eliminate the possibility of 
inverted hierarchy. We should note that uncertainties 
in the z^-process cross sections still remain. We expect, 
however, that they are largely canceled out when we take 
the ^Li/^^B ratio. Since ^Li and ^^B are mainly produced 
through the z^-process from ^He, the dependence of their 
yields on the z^-process reaction rates is similar. In addi- 
tion, the dependence of neutral-current reaction rates on 
the neutrino temperature is not so different from that of 
the corresponding charged-current reactions. Data anal- 



ysis of SN 1987A and future observations of SN relic 
neutrinos may provide additional information on the 
Ue spectrum. The effect of neutrino oscillations on the 
analyzed Pc signal should be taken into account, and the 
evaluation of the Dc spectrum will lead to a more precise 
evaluation of ^Li/^^B ratio in SNe. 

Recent observational efforts to obtain Li and B abun- 
dances in stars which may have formed in regions directly 
affected by prior generations of massive stars and their 
subsequent SNe (e.g., [l^), may have detected the sig- 
nature of the i/-process in ^^B-enriched stars j^^. The 
combination of SN nucleosynthesis theory and observa- 
tions of light elements may ultimately provide powerful 
constraints on mass hierarchy and the mixing angle ^13. 

In summary, we investigated light element synthesis 
in SNe through the i^-process including the change of 
neutrino spectra due to neutrino oscillations. In the case 
of adiabatic 13-mixing resonance and a normal hierarchy, 
the ^Li yield increases by about a factor 1.9 compared 
to the case without oscillations. This increase may be 
accessible to high resolution spectroscopic studies of stars 
in young, star-forming regions. The ^Li yield in other 
cases and the ^^B yield are scarcely affected by neutrino 
oscillations. The adiabaticity of the 13-mixing resonance 
and the mass hierarchy affect robust determinations of 
^Li/"B ratios in SNe. 
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